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ABSTRACT 


Flow  in  the  entrance  region  of  a  parallel  plate  channel  has  been 
theoretically  investigated  by  others  with  little  or  no  experimental  ver¬ 
ification.  Presented  here  is  a  new  method  for  the  analysis  which  in¬ 
volves  a  linearization  of  the  inertia  terms  and  the  introduction  of  a 
stretched  axial  coordinate.  Solutions  are  obtained  for  the  velocity 
distribution  in  an  incompressible  fluid  which  is  valid  from  the  duct 
inlet  to  the  fully  developed  region.  An  experimental  parallel  plate 
channel  was  operated  using  water  as  the  working  fluid.  Velocity  pro¬ 
files,  on  the  basis  of  the  stretched  axial  coordinate,  were  established 
with  the  use  of  a  laser  doppler  velocimeter.  These  velocity  profiles 
are  compared  with  the  theoretical  profiles.  A  significant  variation  be¬ 
tween  the  experimental  and  theoretical  velocity  is  seen  near  the  inlet 
of  the  channel. 
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SECTION  I 
INTRODUCTION 


A  significant  portion  of  the  flow  variations,  particularly  with  short 
length  channels,  is  experienced  in  the  entrance  region.  For  a  constant 
property  fluid,  these  variations  exist  until  ultimately  a  velocity  profile 
is  attained  which  does  not  vary  with  axial  position.  Because  of  the 
influence  of  viscosity,  the  initial  profile  at  the  channel  entrance  will  be 
modified  with  respect  to  axial  position  until  this  fully  developed  profile 
is  achieved. 

A  number  of  theoretical  investigations  of  flow  development  in  the 
entrance  region  of  channels  have  been  made,  but  the  experimental  data 
available  are  meager.  A  summary  of  the  various  analytical  approaches 
to  the  entrance  region  problem  is  presented  by  Sparrow,  Lin,  and 
Lundgren  (Ref.  1).  These  authors  developed  a  new  method  of  flow  anal¬ 
ysis  in  the  entrance  region  and  compared  their  results  with  the  small 
amount  of  experimental  data  which  was  available.  The  only  experimen¬ 
tal  data  cited  were  for  flow  in  a  circular  tube.  The  linearized  analysis 
developed  by  Sparrow  et  al.  gave  better  agreement  with  the  available 
experimental  data  than  the  other  analyses  that  were  cited. 

Recent  experimental  measurements  of  velocity  profiles  in  the  en¬ 
trance  region  of  a  square  duct  have  been  made  by  Goldstein  and  Kreid 
(Ref.  2)  using  a  laser -doppler  flowmeter.  The  experimental  results 
were  compared  with  the  theoretically  calculated  profiles  of  Han  (Ref.  3). 
The  analysis  of  Han  is  based  on  a  Langhaar  (Ref.  4)  method  of  lineari¬ 
zation.  However,  it  has  been  shown  (Ref.  1)  that  the  Langhaar  method 
predicts  a  more  rapid  development  of  the  entrance  region  flow  than  the 
rate  predicted  by  the  linearization  method  when  compared  with  the 
experimentally  measured  rate  of  development.  The  experimental  data 
(Ref.  2)  agreed  well  with  the  theoretical  calculations  in  the  fully  de¬ 
veloped  portion  of  the  channel,  but  variations  existed  between  the  the¬ 
oretical  and  experimental  values  as  one  approached  the  entrance  region. 
The  variations  between  the  theoretical  and  experimental  values  increased 
with  decreasing  distance  to  the  entrance  region. 

In  the  present  investigation,  measurements  of  velocity  profiles  in  the 
entrance  region  of  a  rectangular  channel  were  made  with  a  helium-neon 
gas  laser  anemometer.  The  channel  had  an  aspect  ratio  of  7.  4  to  1.  A 
comparison  was  made  between  the  experimental  measurements  and  pro¬ 
files  calculated  using  the  linearization  method  of  Sparrow,  Lin,  and 
Lundgren  described  in  Ref.  1.  The  principal  objectives  were  to  deter¬ 
mine  the  validity  of  the  analytical  model  described  in  Ref.  1  and  to  .obtain 
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measured  velocity  profiles  in  the  entrance  region  of  a  rectangular 
channel  using  a  laser  anemometer. 


SECTION  II 
ANALYSIS 


In  analyzing  entrance  region  flows,  two  levels  of  approximation 
are  generally  involved,  and  the  distinction  between  them  is  not  always 
made  clear.  The  flow  in  the  entrance  region  is  described  exactly  by 
the  complete  Navier- Stokes  equations  (for  continuum  flow).  The 
Navier-Stokes  equations  are  elliptic,  and  this  means  that  boundary 
conditions  must  be  specified  at  each  point  around  a  closed  contour. 

Thus,  in  addition  to  prescribing  a  profile  at  the  entrance  to  the  channel 
at  x  =  0,  the  profile  must  also  be  specified  at  x  =  ®.  Since  the  boundary- 
layer  equations  are  parabolic,  the  equations  may  be  integrated  in  a  for¬ 
ward  marching  manner  by  starting  with  a  prescribed  initial  profile  at 
x  =  0.  Thus,  the  first  level  of  approximation  involves  the  use  of  the 
bound  ary -layer  form  of  the  equations  of  motion  instead  of  the  full 
Navier-Stokes  equations,  thereby  changing  the  system  from  elliptic  to 
parabolic.  Any  analysis  of  entrance  region  flows  based  on  the  boundary- 
layer  equations  will  be  invalid  at  the  channel  entrance  because  of  the 
neglect  of  the  term  9  2u/  9x2,  which  is  not  small  at  the  entrance. 

The  second  level  of  approximation  is  involved  in  the  particular 
technique  used  to  solve  the  boundary -layer  equations.  Several  approxi¬ 
mate  methods  have  been  used  to  solve  the  boundary -layer  equations;  how¬ 
ever,  in  the  present  analysis  the  approximate  linearization  method  (Ref.  1) 
was  used.  This  method  was  chosen  because  it  appeared  to  give  the  best 
agreement  with  the  available  experimental  data.  The  analysis  is  per¬ 
formed  for  a  parallel  plate  channel  geometry. 

The  geometry  of  the  parallel  plate  channel  is  shown  in  Fig.  1 
(Appendix  I).  The  boundary -layer  form  of  the  equation  of  motion  for  the 
entrance  region  flow  of  a  constant  property,  incompressible  fluid  may  be 
written  as 


9u  du  1  <?p  dJu 

U  <9x  +  '  <9y  p  <?x  V  dy3  (1) 

Because  of  the  nonlinearity  of  Eq.  (1),  an  exact  solution  is  not  possible. 
Following  the  linearization  technique  described  in  Ref.  1,  Eq.  (1)  is 
linearized  as; 
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f(x)U  =  A(x)  -  v 


(2) 


where  e  (x)  is  a  stretching  factor  and  U  is  the  mean  velocity  defined  as 


U  =  ~  /udy 
h  o 


(3) 


The  term  A(x)  contains  the  pressure  gradient  and  the  residual  of  the 
inertia  terms  and  may  be  evaluated  by  integrating  Eq.  (2)  over  the 
channel  half -width  with  the  result 


AW 


Vv_  <9ul 

=  “Lh 


(4) 


Combining  Eqs.  (2)  and  (4)  gives 


/  \TT  (  5*11 

,WU  T,  "  "  |  V 


B-  C ! 


(5) 


A  stretched  coordinate  x*  may  be  introduced  by  the  definition 


dx  =  fdx* 


which  transforms  the  left  side  of  Eq.  (5)  to  a  form  with  a  constant 
coefficient 


cWU 


u 


du 

dx* 


Dimensionless  quantities  may  be  defined  as  follows: 


The  dimensionless  form  of  the  equation  of  motion  becomes 


da  _  d2a 
<5/3  *  drj1 


(6) 


(7) 


(8) 
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with  boundary  conditions 


a(0,i?)  =  a0  (j; );  a(j8,l)  =  0;  ^-(/3;0)  =  0 


(9) 


For  a  uniform  initial  velocity  profile,  Q’0(r])  =  1.  The  solution  to 
Eq.  (7),  subject  to  the  boundary  conditions  of  Eq.  (8)  and  o^fa)  =  1,  is 
given  as  (Ref.  1): 


a 


3 

T 


(1  -  r,1) 


{c°s  Yn V 
COS  Yn 


(-yJ  n/3“) 


(10) 


where  7n  are  eigenvalues  which  are  the  roots  of  the  transcendental 
equation 


Yn  =  tan  yn  (n) 

The  first  98  eigenvalues  of  Eq.  (10)  were  evaluated  (Ref.  5)  and  are 
tabulated  in  Table  I  (Appendix  II). 

To  complete  the  analysis,  the  stretching  factor  e  (x*)  must  be  de¬ 
termined.  The  procedure  developed  in  Ref.  1  is  to  integrate  the  equa¬ 
tion  of  motion  across  the  half -channel  and  then  to  integrate  the  first 
moment  of  the  equation  motion,  i.  e. ,  the  integral  of  the  velocity,  u, 
times  the  equation  of  motion.  Equating  the  pressure  gradient  obtained 
from  the  two  equations  then  results  in  a  solution  for  e  of  the  form 


It  should  be  observed  that  the  evaluation  of  e(|3*)  as  given  by 
Eq.  (12)  is  not  unique  since  other  moments  of  the  equation  of  motion 
could  be  used  as  well.  There  appears  to  be  no  compelling  physical 
reason  for  using  moments  higher  than  the  zeroth  and  first,  and  conse¬ 
quently  these  two  moments  were  used  to  eliminate  the  pressure  gradient 
from  the  equation  of  motion,  resulting  in  a  solution  for  e(/3*). 

The  integrals  appearing  in  Eq.  (12)  were  evaluated  numerically 
using  an  IBM  360/50  digital  computer.  From  the  calculated  values  of 
€(£*),  the  relation  between  the  physical  and  stretched  coordinates  may 
be  obtained  from  the  definition  of  e  as 
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T  =  Re<^*  eifi*  (13) 

Thus,  Eqs.  (9),  (11),  and  (12)  constitute  the  solution  for  the  en¬ 
trance  region  velocity  profile  based  on  the  linearization  model  described 
by  Sparrow  et  al.  (Ref.  1). 

Figure  2  is  a  plot  of  (3  versus  obtained  from  the  solution  of 
Eq.  (13).  Figure  3  is  a  plot  of  a  and  9  or/ 9  j3*  versus  ,6*  evaluated  at  the 
centerline  obtained  from  Eq.  (9).  These  two  quantities  will  be  useful 
in  performing  an  error  analysis  which  will  be  described  later. 


SECTION  III 

DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 


To  determine  the  validity  of  the  above  theoretical  analysis,  an  ex¬ 
perimental  channel  flow  facility  was  designed  and  built  which  would  per¬ 
mit  measurement -of- entrance  region  velocity  profiles  using  a  laser 
anemometer.  The  flow  system  shown  schematically  in  Fig.  4  consists 
of  a  constant  head  reservoir  which  feeds  a  plenum  chamber.  From  the 
plenum  chamber  the  flow  enters  a  2-m-long  developing  section  which 
allows  the  flow  to  become  fully  developed  before  the  initial  conditioning 
section  is  entered.  The  developing  section  was  added  because  of  the 
extreme  difficulty  encountered  in  obtaining  a  constant  initial  velocity 
profile  entering  the  test  section.  Since  the  velocity  profile  existing 
from  this  section  is  the  well-known  Poiseuille  profile,  various  schemes 
for  obtaining  a  flat  profile  at  the  entrance  to  the  test  section  could  be 
studied  with  a  minimum  of  difficulty.  Beyond  the  developing  section  the 
flow  passes  through  a  screen  diffuser  with  a  solidity  ratio  of  0.  162, 
i.  e. ,  1  -  Aq/A-j.  when  AQ  is  the  open  area  of  the  screen  and  Af.  is  the 
total  area  of  the  screen,  and  into  the  test  section  which  consists  of  two 
rectangular  glass  channels,  one  inside  the  other.  The  dimensions  are 
14  by  3.  8  cm  for  the  outer  channel  and  14  by  1.  9  cm  for  the  inner  chan¬ 
nel.  The  purpose  of  the  inner  channel  is  to  allow  only  the  flat  central 
portion  of  the  channel  profile  to  enter  the  test  section,  thus  resulting 
in  a  well-defined  flat  initial  profile.  The  length  of  the  test  section  is 
76  cm.  The  leading  edges  of  the  inner  channel  are  slightly  rounded  to 
minimize  flow  separation.  The  data  presented  were  for  the  flow  within 
the  inside  channel.  The  flow  from  the  test  section  then  enters  a  dump 
tank  where  it  is  returned  to  the  reservoir  by  the  secondary  loop  shown 
in  Fig.  4.  It  was  found  necessary  to  vent  the  flow  system  to  remove 
fluctuations  caused  by  air  bubbles  trapped  during  the  filling  process. 
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The  laser  anemometer  system,  shown  schematically  in  Fig.  1,  has 
a  traverse  capability  across  the  channel  and  down  the  length  of  the  chan¬ 
nel  so  that  velocity  profiles  can  be  obtained  at  different  stages  of  flow 
development.  The  anemometer  incorporates  the  doppler  shift  effect 
with  optical  homodyne  detection  similar  to  that  used  by  Yeh  and  Cummins 
{Ref.  6). 


SECTION  IV 
INITIAL  CONDITION 


Considerable  effort  was  expended  to  obtain  a  flat  velocity  profile  at 
the  beginning  of  the  knife  edges  of  the  test  channel.  Figure  5  is  a  super¬ 
imposed  plot  of  the  velocity  profile  immediately  upstream  and  downstream 
of  the  straightening  screen.  The  upstream  profile  is  parabolic,  as  would 
be  expected  since  the  flow  has  passed  through  the  2-m  developing  section. 
The  downstream  profile  exhibits  a  marked  concavity  on  the  centerline.  It 
was  found  that  a  very  minute  deflection  in  the  center  of  the  straightening 
screen  was  the  apparent  reason  for  the  inflection.  The  screen  holder  was 
redesigned  so  that  the  screen  could  be  put  under  considerable  tension  to 
avoid  the  deflection.  The  resulting  profile  downstream  of  the  new  screen 
arrangement  is  shown  in  Fig.  6.  This  is  a  typical  initial  profile  obtained 
for  all  the  experiments  that  were  conducted.  The  location  of  the  knife 
edges  relative  to  the  profile  is  also  shown  in  Fig.  6. 


SECTION  V 

ANALYSIS  OF  THE  EXPERIMENTAL  DATA 


Since  the  doppler  frequency  is  directly  proportional  to  the  longitu¬ 
dinal  component  of  velocity,  for  the  optical  configuration  used  in  this 
study,  normalized  velocity  profiles  were  obtained  by  numerical  integra¬ 
tion  of  the  doppler  frequency  shift  profiles.  The  normalized  profiles 
across  the  channel  were  obtained  by  dividing  the  locally  measured  fre¬ 
quency  by  the  average  frequency  over  the  channel  width  and  normalized 
to  the  center  velocity  at  the  entrance  of  the  channel;  hence,  the  center, 
fully  developed  Poiseuille  velocity  becomes  1.5.  Figure  7  is  a  plot  of 
several  normalized  velocity  profiles  obtained  in  this  manner  using  j3*  as 
a  parameter.  By  using  j8*  as  a  parameter,  one  may  make  all  measure¬ 
ments  at  one  downstream  position,  x,  and  vary  the  flow  rate  {Reynolds 
number)  to  produce  the  equivalent  developing  profile.  The  Reynolds 
number  used  to  calculate  j8*  was  obtained  from  flow  rate  measurements 
determined  by  volumetric  means.  It  is  apparent  that  the  measured  and 
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calculated  profiles  are  -in  close  agreement  near  the  fully  developed 
region  (large  values  of  /3*).  Near  the  entrance,  however,  the  meas¬ 
ured  centerline  velocities  are  consistently  lower  than  the  calculated 
values.  To  check  the  value  of  ,3*  obtained  from  the  volumetric  meas¬ 
urements,  the  average  doppler  frequency  was  converted  to  the  flow  rate 
and  compared  with  the  volumetric  flow  rate.  The  doppler  frequency  is 
given  by 


2n v  .  e  e 

-  sin  —  cos  — 

A0  2  2 


(14) 


where  fj)  is  the  frequency  measured,  V  is  the  velocity,  n  is  the  index 
of  refraction,  and  6  is  the  angle  between  scattered  beam  and  reference 
beam.  The  laser  wavelength  is  A0.  The  mass  flow  rate  W  is  given  by 


W  =  <V>  A 


(15) 


with  A  being  the  cross-sectional  area  of  the  test  section  and  where  <V> 
is  the  average  velocity  or  bulk  velocity. 

Since  the  average  velocity  is  directly  proportional  to  the  average 
doppler  shifted  frequency,  W  may  be  found.  The  average  frequency  was 
established  from  a  computer  program  which  provided  an  analytic  equa¬ 
tion  for  each  velocity  profile  (frequency  profile)  and  performed  the 
following  integration. 


<f>  =  -1  fffdA 
A 


With  the  use  of  the  average  frequency  in  Eq.  (14)  and  substituting  into 
Eq.  (15),  we  have 


<fD>  A0  A 

in  (f)  cos 


2n  sin 


(16) 


The  angle  9  was  found  with  the  use  of  Snell's  law  (also  see  Fig.  1) 


sin  0!  =  na  sin  d3 


(17) 


where  nj  and  n2  were  the  index  of  refraction  of  the  water  and  air,  re¬ 
spectively.  Since  the  data  were  taken  over  a  period  of  several  days,  a 
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new  0 2  was  found  each  day  to  ensure  high  accuracies.  Then  8 1  was 
calculated  with  n  =  1  and  n£  =  1.  35  such  that 


6X  =  sin-1 


(18) 


Table  II  contains  the  flow  rate  obtained  in  this  manner  along  with 
those  obtained  from  the  volumetric  measurements.  Except  for  run  1,  it 
is  seen  that  the  laser  anemometer  values  are  consistently  lower  than 
the  volumetric  measured  values.  Hence,  the  values  of  0*  calculated 
from  Reynolds  numbers  based  on  the  anemometer  measurements  would 
then  be  larger  than  those  from  the  volumetric  measurements  by  an 
average  of  2.  8  percent.  Since  the  normalized  centerline  velocity  in¬ 
creases  with  increasing  j3*,  profile  comparisons  based  on  the  anemom¬ 
eter  /3*'s  would  show  an  even  worse  correlation  near  the  entrance.  The 
shape  of  the  measured  profile  indicates  that  the  disagreement  might  be 
explained  in  terms  of  an  error  in  the  flow  rate  measurement,  which 
would  in  turn  produce  an  error  in  |3*. 

To  check  this  possibility,  an  error  analysis  of  the  centerline  veloc¬ 
ity  determination  attributable  to  an  error  in  the  measured  flow  rate  was 
performed.  The  error  in  centerline  velocity  owing  to  an  error  in  /3*  is 
given  by 


From  the  definition  of  0*,  Eq.  (7),  the  error  in  the  calculated  /3*  pro¬ 
duced  by  an  error  in  measured  flow  rate  is  given  by 

ARe  „  Aw 

A/3*  0-/3*-^ - 0*  -j-  (20) 


providing  the  increment  is  small  compared  to  the  absolute  magnitude. 
Thus, 


(Aa)ci  = 


(21) 


The  derivative  was  found  by  plotting  j3*  versus  the  experimental  center- 
line  velocity  and  then  performing  a  graphical  differentiation  of  this  plot. 
The  results  are  reproduced  in  Fig.  3. 
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The  percentages  of  error  between  the  theoretical  and  experimental 
data  of  the  /3*1  s  associated  with  the  normalized  velocity  profiles  shown 
in  Fig.  7  are  tabulated  in  Table  III.  In  addition,  the  table  includes  the 
error  between  the  anemometer  and  volumetric  flow  rates,  assuming  the 
anemometer  measurements  were  correct.  These  experimental  flow 
rate  errors  were  used  in  Eq.  (21)  to  determine  the  approximate  error 
in  the  centerline  velocity.  The  results  are  tabulated  in  the  last  column 
of  Table  III.  It  is  seen  that  this  error  is  negligible  compared  with  the 
error  between  the  experimental  and  theoretical  data. 

Several  sources  of  error  are  attributable  to  the  laser  anemometer, 
but  these  errors  may  be  assumed  to  be  quite  small.  These  errors  are 
inherent  in  determining  the  angle  02  (°r  the  ability  to  measure  the  three 
sides  of  a  triangle  accurately),  the  ability  of  determining  the  mean  veloc¬ 
ity  of  a  Gaussian  velocity  distribution,  and  the  ability  to  measure  fre¬ 
quency  accurately.  These  errors  were  found  to  be  less  than  0.  5  percent; 
hence,  they  contributed  little  or  no  error  in  ,3*.  The  location  of  the  posi¬ 
tion  x*  was  found  within  0.  0001  in. 

At  this  time,  it  may  be  assumed  that  the  gradual  increase  in  varia¬ 
tion  that  is  found  by  approaching  the  channel  entrance  is  caused  by  insuf¬ 
ficiencies  in  the  theoretical  model  on  which  the  analysis  is  performed. 


SECTION  VI 
CONCLUSIONS 


The  velocity  profiles  in  the  entrance  region  of  a  rectangular  flow 
channel  were  measured  with  a  laser  anemometer.  The  measurements 
were  compared  with  profiles  calculated  by  a  method  developed  by 
Sparrow.  Satisfactory  agreement  was  found  near  the  fully  developed 
region,  but  deviations  were  noted  close  to  the  entrance  which  cannot 
be  accounted  for  by  experimental  error.  It  is  concluded  from  the  exper¬ 
imental  data  that  the  theory  used  accurately  predicts  entrance  region 
profiles  for  distances  from  the  channel  entrance  such  that  x/h  >  (0.025/Re). 
For  distances  smaller  than  this,  the  theory  predicts  a  more  rapid  de¬ 
velopment  of  the  flow  than  is  measured  experimentally.  This  is  thought 
to  be  an  inadequacy  in  the  analytical  model  very  near  the  entrance.  Since 
the  model  used  in  the  analysis  is  of  the  boundary -layer  type,  a  discrep¬ 
ancy  between  theory  and  experiment  would  be  expected  very  near  the 
entrance  because  of  the  neglect  of  the  9^u/9x^  term  in  the  boundary- 
layer  equation. 


9 


AEDC-TR-69-186 


REFERENCES 


1.  Sparrow,  E.  M.  ,  Lin,  S.  H. ,  and  Lundgren,  T.  S.  "Flow  Develop¬ 

ment  in  the  Hydrodynamic  Entrance  Region  of  Tubes  and  Ducts, 11 
Physics  of  Fluids.  Vol.  7,  No.  3,  March  1964,  pp.  338-347. 

2.  Goldstein,  R.  J.  and  Kreid,  D.  K.  "Measurement  of  Laminar  Flow 

Development  in  Square  Duct  Using  Laser-Doppler  Flowmeter,  " 
Journal  of  Applied  Mechanics.  December  1967,  pp.  813-818. 

3.  Han,  L.  S.  "Hydrodynamic  Entrance  Lengths  for  Incompressible 

Laminar  Flow  in  Rectangular  Ducts,  "  Journal  of  Applied 
Mechanics.  September  I960,  pp.  403-409. 

4.  Langhaar,  H.  L.  "Steady  Flow  in  the  Transition  Length  of  a  Straight 

Tube,"  Journal  of  Applied  Mechanics.  June  1942,  pp.  55-58. 

5.  Snyder,  W.  T.  and  Eraslan,  A.  H.  "The  Influence  of  Hall  Effect 

and  Initial  Velocity  Profile  on  MHD  Flow  in  the  Entrance 
Region  of  a  Parallel  Plate  Channel.  "  AEDC-TR-67-79 
(AD650820),  April  1967. 

6.  Yeh,  Y.  and  Cummins,  H.  Z.  "Localized  Fluid  Plow  Measurements 

with  He-Ne  Laser  Spectrometer,  "  Applied  Physics  Letters. 

Vol.  4,  No.  10,  May  1964,  pp.  176-178. 


10 


A  EDC-TR-69- 186 


APPENDIXES 

I.  ILLUSTRATIONS 

II.  TABLES 


11 


Diffuser 
Unit  — 


Screen 


AEDC-TR.69.189 


Fig.  4  Schematic  Diagram  of  Flow  System 
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Fig.  5  Velocity  Profiles  of  Upstream  and  Downstream  Sides  of  a  Slightly  Deformed 
Diffuser  Screen 


Normalized  Velocity,  u/u 
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Fig.  7  Normoliied  Velocity  Profiles  versus  Normalized  Position 
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TABLE  I 

EIGENVALUES  OF  TAN  y„  =  yn 


n 

7n 

n 

7n 

n 

Vn 

— 

1 

4.4934127 

33 

105. 23392 

66 

208. 91127 

2 

7. 7252596 

34 

108.  37579 

67 

212. 05293 

3 

10. 904133 

35 

111.  51764 

68 

215. 19459 

4 

14.066202 

36 

114. 65947 

69 

218.  33625 

5 

17.220769 

37 

117.  80130 

70 

221.47790 

6 

20.371314 

3& 

120.  94311 

71 

224. 61956 

7 

23.  519473 

39 

124. 08491 

72 

227. 76121 

8 

26. 666072 

40 

127.  22676 

73 

230.  90286 

9 

29.  811615 

41 

130.  36854 

74 

234. 04451 

10 

32. 956418 

42 

133.  51031 

75 

237.  18616 

11 

36. 100648 

43 

136. 65207 

76 

240. 32781 

12 

39.244457 

44 

139.  79383 

77 

243.46945 

13 

42.387936 

45 

142.  93557 

78 

246.61109 

14 

45. 531176 

46 

146.07731 

79 

249.  75286 

15 

48.674183 

47 

149.  21905 

80 

252.  89450 

16 

51.817019 

48 

152. 36077 

81 

256.03613 

17 

54. 959713 

49 

155. 50250 

82 

259.  17777 

18 

58. 102287 

50 

158. 64422 

83 

262.31941 

19 

61.244761 

51 

161.  78593 

84 

265.46104 

20 

64.  387179 

52 

164. 92764 

85 

268. 60268 

21 

67. 529490 

53 

168.  06935 

86 

271.  74431 

22 

70. 671740 

54 

171.  21104 

87 

274.88594 

23 

73.813931 

55 

174. 35274 

88 

278.  02757 

24 

76. 956076 

56 

177.49453 

89 

281.  16920 

25 

80.098176 

57 

180. 63621 

90 

284. 31083 

26 

83.240237 

58 

183. 77789 

91 

287.45246 

27 

86. 382265 

59 

186.  91958 

92 

290. 59408 

28 

89. 524307 

60 

190. 06126 

93 

293. 73571 

29 

92.666274 

61 

193. 20293 

94 

296. 87734 

30 

95. 808220 

62 

196. 34460 

95 

300. 01897 

31 

98.  950141 

63 

199.48627 

96 

303. 16059 

32 

102. 09204 

64 

65 

202.  62794 

205. 76961 

97 

98 

306. 30222 
309.44384 
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TABLE  II 

COMPARISON  BETWEEN  MASS  FLOW  RATES  OF  ANEMOMETER 
AND  YOLUMETRIC  MEASUREMENT  AT  APPROXIMATELY  85CF 


Run  No. 

Flow  Rate, 
cm3 /sec 
(Anemometer) 

Flow  Rate, 
cm3/sec 
(Volumetric) 

Deviation  r~ ^ »  percent 

1 

16.  0 

15.  7 

-1.  91 

2 

58.  8 

60.  2 

2.  33 

3 

15.  5 

15.  9 

2.  52 

4 

59.  1 

60.  7 

2.  64 

5 

56.  3 

60.  2 

6.48 

,  6 

33.  5 

34.  4 

2.  62 

7 

50.  9 

52.  7 

3.42 

8 

72.  0 

73.0 

1.  37 

9 

64.  5 

66.  5 

3.  01 

10 

119.  5 

127.  0 

5.  91 

„  (V  -  A)  =  2.  84  percent 
i"1 
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TABLE  III 

PERCENT  OF  ERROR  IN  DATA 


Error,  percent  Error,  percent  Error,  percent 


(Experimental  -  Theoretical) 

(Volumetric  -  Anemometer) 

(Experimental  o-cl) 

0.440 

0 

2.52 

~0 

0.0892 

1.  81 

2.62 

0.44 

0.0491 

1.56 

3.42 

0.006 

0. 0295 

7.  39 

1.37 

0.  017 

0.0175 

5.  75 

6.  72 

0.006 
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